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Abstract 

Using e-expansion technique we compute rj/s, where n is the shear viscosity, s is the entropy 
density, of the normal phase of unitary Fermi gas in d = 4 — e dimensions to LO in e. We use kinetic 
theory approach and solve transport equations for medium perturbed by a shear hydrodynamic 
flow. The collision integrals are calculated to e 2 which is LO. The LO result is temperature 
independent with 77/s ~ (0.11/£ 2 )(fo/kB). The d = 3 prediction for 77/s exceeds the h/4nkB bound 
by a factor of about 1.4. 
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I. INTRODUCTION 



Cold dilute gas consisting of two Fermion species (which will be labeled as |, 4) with two 
body interactions characterized by 

no 3 > 1, nrl < 1, 

where n is the number density, a is the s wave scattering length and ro is the effective range 
of the potential, named "unitary Fermi gas" is one of the simplest examples of a strongly 
interacting Fermi system. Nevertheless, the unitary Fermi gas is only a slight idealization 
of dilute neutron matter (a ~ — 18/m, r — 2.6/m in the 1 So channel) that may exist in 
the crust of neutron stars [1, 2]. Also, it has been studied extensively in the experiments on 
trapped cold atomic gases with tunable interactions using the Feshbach resonance technique, 
e.g., [3-7]. In these experiments, they have an amazing control over both the strength of the 
interaction and the composition (i.e., imbalance of concentrations of pairing Fermi particles, 
or polarization) of the system (see [8] for a review of experiments). From the theoretical 
point of view, the unitary Fermi gas lacks any intrinsic scale and is expected to exhibit 
universal properties. At the same time, a theoretical description of unitary Fermi gas is 
difficult due to the apparent absence of a small dimensionless parameter which could be 
used in a perturbative expansion, thus making the usual perturbative techniques unreliable. 

Properties of these systems have been studied in various Monte Carlo simulations (see [9] 
for a compilation of results), and using various theoretical approaches ranging from large- 
N expansion [10] and functional renormalization group [11] to phenomenological density 
functionals [12]. See [13] for a review. But an analytical systematically improvable technique 
is still needed. For example, a Monte Carlo approach is not well-suited to study dynamical 
quantities (such as dynamical response function) and becomes impractical for the polarized 
case due to the sign problem. 

Inspired by the ideas of Nussinov and Nussinov [14], an analytical technique similar to 
the e-expansion in the theory of critical phenomena [15] has been proposed by Nishida and 
Son [16]. The idea is to describe three dimensional unitary Fermi gas using perturbation 
theory in the dimensionality of space, d, around d = 4 and d = 2. The rational for this 
approach is that unitary Fermi gas simplifies to 

• non-interacting Fermi gas in d — 2, 
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• the system of non-interacting spin dimers in d = 4. 

This approach is also applicable to the case of large but finite scattering length (Fermi 
gas near unitarity) [17, 18]. Predictions for d = 3 may be obtained by extrapolating the 
series expansion to e — 1, possibly employing various convergence improvement techniques, 
such as Pade approximation to the Borel transformation of the series, conformal mapping 
technique, etc. [19]. 

This approach has been used to study zero temperature equation of state to the next-to- 
next-to-leading order (NNLO) in e [20, 21]. For the Bertsch parameter, £, use of expansions 
around d — 4 and d = 2 with the use of Pade approximants has given £ = 0.36 ±0.02 [21]. A 
recent fixed-node diffusion Monte Carlo simulation has yielded £ < 0.383(1) [22] and recent 
lattice calculations extrapolated to zero temperature have produced £ = 0.37(5) [23] and 
e = 0.39681™™ [24]. 

To NLO in e the expansion around d — 4 has been used to calculate the fermion quasi- 
particle spectrum [16, 18] and good agreement with Monte Carlo results was found [25]. See 
[26] for a review. Also, atom-dimer and dimer-dimer scattering in vacuum has been studied 
[27]. The low energy density response function has been calculated [28]. The polarized 
case has been studied and a polarization induced superfluid-normal phase transition has 
been found [18, 29]. Also, effective Lagrangian functional (Landau-Ginzburg-like functional) 
for this system has been derived and used to study the structure of a superfluid vortex, 
and of the superfluid-normal interface in the polarized case, all to NLO [30]. The finite 
temperature properties of the unitary Fermi gas have been studied to NLO. In particular, 
the critical temperature of the superfluid-normal phase transition has been calculated and 
good quantitative agreement with numerical results has been found [31]. 

In all of the above, the NLO corrections have been found to be reasonably small, especially 
in the non-polarized case. All these results suggest that already at the next-to-leading order 
(NLO) in e the expansion around d — 4 may be a useful tool for the unitary Fermi gas 
observables, and that working to a higher order in e combined with employing results of 
expansion around d = 2 and convergence improvement techniques is expected to produce 
improvements in accuracy. 

So, as another step in the investigation of the e-expansion usefulness here we compute 
shear viscosity of the normal phase of unitary Fermi gas in d = 4 — e spatial dimensions to 
LO in e. Shear viscosity, 77, is a characteristic of the fluid's internal friction. Determination 
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of this particular transport coefficient has received a lot of attention in recent years, e.g., 
[32]. See [33] for a review. This is, in part, due to the fact that in the strong coupling limit 
of a large class of systems with gauge interactions that can be analyzed from the AdS/CFT 
correspondence rj/s — h/Anks, where s is the entropy density, ks is the Boltzmann constant. 
And a lower bound 



for any substance has been proposed [34]. Both the quark gluon plasma produced in the 
relativistic heavy ion collisions and the unitary Fermi gas have been shown to have [35-38] 



and are considered as the perfect fluid candidates. In this work (1) will be treated as a 
benchmark allowing one to determine usefulness of the e-expansion technique. 

In Section II shear viscosity calculation is set up in the e-expansion framework. In Section 
III calculation of the normal phase collision integrals in d = 4 — e to LO in e is presented. In 
Section IV the linearized transport equation solutions for the shear hydrodynamic flow are 
presented. Using the previously calculated result for the normal phase pressure LO result for 
i]/s is presented which is then compared with the rj/s predictions by some other methods. 

II. LINEARIZED TRANSPORT EQUATIONS AND EXPRESSION OF THE 
SHEAR VISCOSITY 

Due to universality of the unitary Fermi gas, we choose to work with local four-Fermi 
interaction with a coupling constant Co tuned to reproduce infinite scattering length. The 
Lagrangian density of the unitary Fermi gas is (here and below K = k B = 1) 



where ip = (^,^) T , is spin-1/2 fermion field, and = ipj ipi, with % = { |, | } are particle 
density operators, /x is the particle number chemical potential which in this calculation is 
assumed to be the same for both pairing species. The theory is defined in d = 4 — e spatial 
dimensions. After Hubbard- Stratonovich transformation, the Lagrangian density (3) takes 



r)/s < 0.5 



h 

k B 



(2) 




(3) 
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the form 



-— 0*0 + # t <7+#0 + * f (7_*0*, (4) 

Co 

where is an auxiliary di-fermion field, \1/ is a two-component Nambu-Gor'kov field, 
^ = (^,^) T , a± = \{<J\ ± 202)) with o"i,2,3 being the Pauli matrices. In this calculation 
we will concentrate on the unitary regime, a = 00, which in the dimensional regularization 
corresponds to l/c = [16]. Once fluctuations of are identified as propagating bosonic 
degrees of freedom, the system may be viewed as spin up and spin down fermions inter- 
acting by exchanging spin mass 2m bosons. The standard perturbative approach applies 
[16, 18]. Extension to the finite temperature-real time perturbation theory (e.g., [39]) is 
straightforward. The fermion-boson coupling for the LO calculations is [16] 

9 = [ ^^—. (5) 
m 

Transport coefficients can, in principle, be determined using Kubo relations of linear 
response theory. However, it has been argued that using a perturbation theory may require 
re-summation. The equivalent correct procedure is to solve the appropriately perturbed 
Boltzmann equation in which the collision integral is calculated to a finite order in the 
perturbation theory [40]. This is the approach used in this work. 

Momentum flux tensor of a weakly non-ideal non-relativistic fluid is [39] 

n Q(3 = P5 Q/ 3 + pV a V,g - (r' ap , a, /3 = l..d, (6) 

where P is pressure, V(x) is the hydrodynamic flow velocity field, a' a/3 is the viscous stress 
tensor. For a static shear flow V • V = 

1 (dNn <9V n 



^ = 2^,V„, = -^ + ^j,V^0, (7) 

where 77 is the (static) shear viscosity. 

We have three interacting liquids: spin up and spin down fermion species interact with 
each other by exchanging bosons. In the unpolarized case we are considering here there are 
two distribution functions = = n F , n B . For a static shear flow the transport equations 
become 



— = -n F (\ - n F ) \ PaPfs — \-^f = C F [n F ,n B \, 



(8) 



^ = -n° (1 + n%) (p aP , - ^ = C B [ n „n B ], (9) 



where % are the Boson and Fermion distribution functions, respectively, and C F , C b 
are the collision integrals, and 

n% = (exp(E F (p)/T) + l)" 1 , n% = (exp(E B (p)/T) - l)" 1 , 
Ef(p) = e p -n, E B (p) = e p /2 - 2/x, e p = p 2 /2m, (10) 

are the equilibrium Fermion and Boson distribution functions, respectively [39]. One seeks 
solutions as 

n F {p) = rip- n F (l - n F ) PaPp -~rP ] 



d I mT 2 



( \ 0/1 , \ / <W 2 \ Xfi( e p) /,,\ 

«B(P) =»B- WbI 1 + U b) [PaPf) ~-fP ) ■ (11) 

Subsituting (11) into transport equations (8), (9) and retaining terms linear in xf, Xb one 
obtains two coupled linea integral equations for xf, Xb- Once xf, Xb have been determined 
i] is given by 



III. LO COLLISION INTEGRALS 

It is known that in the finite temperature-real time formalism the equation of motion 
for the two point correlation function G h (t 1 , x 1; t 2 , x 2 ) = ±i(^£ fa, X2)V'a(ii, x i))> where ± 
refers to Fermions and Bosons, respectively, in the quasiclassical, or Thomas-Fermi, regime 
reduces to the transport equation for the distribution functions [39]. The LO collision 
integrals can then be read off and are given by 

C B [n F ,n B ] = ?S B + (1 + n B ) - iS^n B | fco=ek/ 2-2 M , 

C F [n F ,n B } = «S^ + (1 - n F ) +iE+~n F | feo=£k _ jtii (13) 
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where iS^(/c , k; rip, n B ), i£^(/c , k; n F , n B ) are the boson and fermion self energies [39]. 
In Eq. (13) ns may depend on time as well as on the phase space coordinates, x and k, 
which are assumed to be slow varying on the microscopic scales. 

The LO contributions to the collision integral are 0(e 2 ). This is as expected since the 
LO contributions to fermion-fermion and fermion-boson scattering cross-sections are all 
0(e 2 ). We have checked that the 0(e) contributions to both boson and fermion self-energies 
iS^~(/c , k), iS^ + (/c , k) and iS^ _ (/c , k), iS^ + (/c , k) vanish for k = e k /2 — 2/i, k = e k — /i, 
respectively. 

Feynman diagrams relevant for the calculation of the collision integrals in the normal 
phase are shown in Fig. (1). For the rest of this section we will be distinguishing f and J, 
fermion species. This is just for convenience. At 0(e 2 ) the fermion collision integral is 

C t4 = Cj 4 + C? 4 , (14) 

where 

c t = / m p an q x 

x [n B (q + 2k) (1 + n B (p + 2q)) n t (p + q)(l - n t (k)) - 

- (1 + n B (q + 2k)) n B (p + 2q)(l - n t (p + q))n t (k)}, (15) 

and 

a 4 (2m) d f 

cf = J de q <m p <m q x 

x Kb + q)Mq - p)(! - n ^( 2 9 + fe))(i - ^t(fc)) - 

- (l-n ; (p + q))(l-n t (g-p)K(2g + fc))n t (fc)], (16) 

where \p\ = \q\, q = q + k. Expressions for Cj and are given by formulas (15), (16) 
with | subscripts replaced by | subscripts. C 1 comes from the self energy correction to the 
fermion propagator, while C 2 comes the self energy correction to the boson propagator in 
the fermion self energy diagrams shown in Fig. (1). The 0(e 2 ) contribution to the boson 
collision integral comes from the self energy corrections to the fermion propagators in the 
boson self energy diagram shown in Fig. (1). It is given by 

c B = 9 J d £q <m p <m q x 

x ej- 4 [(l - H (k/2 - q)) H (p + q)(l + n B (p + 2q))n B (k)) - 

- ni (k/2 -q)(l- ni (p + q))n B (p + 2q)(l + n B (k))} + { t<-4 }, (17) 
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FIG. 1: Feynman diagrams relevant for the normal phase calculation the collision integrals. Thin 
solid lines depict propagating fermions, dashed line is the boson propagator. The boson self-energy 
diagram with the boson line in the lower part of the fermion loop is not shown. The — h diagrams 
are completely analogous. The rest of the 0(e 2 ) contributions to the self energies produce zero 
contributions to the collision integrals. 

where |p| = 2|qr|, q = q + k/2. In the LO expressions above d = 4 and, so, 

dVt p — d(f) du sin 2 w / d8 sin6 
Jo Jo Jo 

is the four dimensional angular integration. Also, we note that collision integrals have the 
standard gain-loss structure [39]. Evaluated on the equilibrium distribution functions (10) 
(15), (16), (17) vanish, as they should. 

IV. SOLUTION OF LINEARIZED TRANSPORT EQUATIONS, RESULT AND 
CONCLUSION 

Substituting (11) into (8), (9) two linear integral equations for XF^k), Xs( e k) have been 
obtained. The expressions for the linearized integral operators are somewhat bulky and will 
not be shown here. 

First, we have checked the integrals for several suspected collinear divergences coming 
from the Bose distribution functions when p, q and k from (15), (16), (17) are (anti)parallel. 
We have observed that in d — 4 all the suspected divergences produce finite contributions. 
This has been verified a posteriori numerically using the solutions obtained. 

The linearized equations have been solved numerically for the hydrodynamic flow with 
V34 7^ (cf. [41]) by denning the unknown functions on a discretized finite domain. Thus 




FIG. 2: Solutions to the linearized transport equations obtained by the method described in text; 
Xb is the solid line, xf is the dashed line. 

the problem has been reduced to the system of linear algebraic equations. However, due 
to non-locality of the kernels (15), (16), (17) the resulting system of linear equations is 
not closed, i.e., for a fixed domain of the argument < e^/T < A the integral operators 
will always depend on values of XF^k), Xs( e k) f° r £k/T > A. An ad hoc solution used 
here is to simply set XF^k), Xs( e k) = for > A. The seven fold integrations have 
been done numerically by the Mathematica software using Adaptive Monte Carlo method 
in the NIntegrate function. We have checked for convergence of our solutions with respect 
to both the lattice spacing and the value of the cutoff, A. The solutions obtained have been 
converted back into the "continuous" form using Interpolation function of Mathematica and 
fed back into the linearized transport equations. We have found that for < e^/T < 13 
relevant to the r\ calculation the ratio of the L.H.S. to the R.H.S. was 1.07 ± 0.07 for the 
linearized Boltzmann equation for fermions and 0.97 ± 0.03 for the linearized Boltzmann 
equation for the boson distribution function. From this it was concluded that reasonably 
accurate solutions had been found. The errors may be further reduced by making the lattice 
finer and the cutoff A bigger, but even at this level they do not seem to be the dominant 
source of uncertainty in the result. Functions Xf{^v)i Xs( e k) are plotted in Fig. 2. Note 
that the approach used here differs from the standard variational approach described in the 
literature where one seeks the solution(s) as an expansion in some basis which results in the 
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lower bound for the viscosity [39]. For some recent applications of the variational approach 
see, e.g., [41, 42]. 

The pressure of the normal phase for d = 4 — e has been computed by Nishida to NLO 
[43] 

P = tY— 9 1112 C(3) + 11^(3) ttV 27^ \ 

V 2tt y V 2 4 p 6T 3 T )' 

C p ~ 8.4, n B = 2n~ e. (18) 

The entropy density is s = dP/dT. Evaluating (12) on the linearized transport equations 
solutions X-F^k), XB^k) and using (18) to LO in e we get (recall h = k B = 1) 

n 1 1 1 44 1 

^^(l + <9( £ )) = — -(l + 0( £ )). (19) 

The LO normal phase ratio turns out to be temperature independent. 

Here we make the simplest possible extrapolation for d — 3 by setting e — 1 in (19). We 
see that our prediction does not contradict the experimentally obtained bound (2) and is 
consistent with the conjectured bound (1). A recent first principles calculation of the normal 
phase rj/s used Kubo linear response approach with the stress tensor correlator computed 
on the lattice [44]. The d — 3 prediction of this work, (19) with e — 1, appears to coincide 
with the (rf/s) min ~ 0.11 quoted in [44]. We believe this to be a fluke. Other theoretical 
approaches yielded r//s ratios exceeding the bound (1) by, e.g., ~ 7 [45] and 4.7 [46]. 

The main conclusion of this work is that e-expansion is a versatile useful tool for the 
unitary Fermi gas already at LO and NLO. Further investigation aimed at high precision 
calculations of various observables is fully warranted. 
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